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SUMMARY 

The ^transonic similarity law for two-dimensional flow derived 
hy von Karman has "been investigated "by an iteration procedure 
similar to that of the Rayleigh-Janzen and Ackeret-Prandtl-Glauert 
methods. The results, which show that the potential can he expressed 
as a power series in a single parameter that depends on Mach number, 
thickness ratio, and ratio of specific heats, are in agreement with 
those, of von Karman. The iteration procedure ha« been applied to 
obtain the second approximation for the flow past a Kaplan section 
in similarity form. The exact solution by Kaplan for the second 
approximation has been examined and found expressible in the same 
similarity form. The exact numerical results to three approximations 
obtained by Kaplan for the Kaplan section and the circular arc have 
been reduced to transonic similarity form. 


INTRODUCTION 

The difficulty of calculating and measuring explicitly fluid 
flow patterns past given bodies has led to attempts at extracting 
information from the equations of motion without actually solving 
them, that Is, in the form of similarity laws. A similarity law 
gives the solution for a whole class of related bodies at related 
conditions if the solution for only one of these bodies at one set 
of conditions is known. A well-known example of a similarity law 
is the Prandtl-Glauert rule, which relates potential flows past thin 
airfoils at low subsonic speeds. 

Similarity rules for two-dimensional steady-state potential • 
flow in the transonic speed range have been derived by von Karman 
in reference 1. The treatment has been extended by Lin, Reissner, 
and Tslen in reference 2 to include unsteady two-dimensional 
transonic potential flows. These derivations utilized small- 
perturbation methods in the physical plane in which two basic 
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simplifying assumptions were made: namely, (1) that the effeot of 
stagnation points, if these exist in the flow field, was negligible 
in the region of interest in the flow field; and (2) that the 
boundary condition of zero velocity component normal to the body 
contour could be satisfied near the body rather than on the body 
itself. The first assumption is also made herein. The validity of 
the second assumption was Investigated at the HA.CA Lewis laboratory 
and is reported -herein. 

A rigorous Justification of the assumption that the exact 
boundary condition may be replaced by another one would appear to 
require the study of the behavior under these various boundary 
conditions of an explicit solution of the equations of motion. An 
explicit solution, however, is Just what the similarity analysis 
purports to do without. A dilemma is thus apparent. That this 
dilemma is nontrivial is illustrated by the fact that the preceding 
boundary-condition assumption has been known to lead to substan- 
tially incorrect results in the case of slender bodies of 
revolution in the subsonic speed range (see, for example, refer- 
ences 3 to 6), although, as is well known, this general assump- 
tion yields correct results in this speed range for two- 
dimensional airfoils. 

The explicit class of solutions that may therefore be needed 
to establish a similarity law does not, however, render superfluous 
the similarity law contained implicitly therein. The expression of 
these solutions in a form that brings out the similarity features 
is both useful and significant, as is evident from the example of 
the Prandtl-Glauert rule. Accordingly, an explicit solution for 
the two-dimensional, continuous, potential flow past thin airfoils 
at high subsonic speeds is herein rigorously derived from the point 
of view of the possibility of expressing this solution in the form 
required by the transonic similarity law (reference 1) . The 
l imi t ing form of the solution is found to be expressible in this 
form as the airfoil thickness ratio approaches zero and the free- 
strearn Mach number approaches unity. Hence, the boundary-condition 
assumption in question, underlying the derivations of references 1 
and 2, is validated for the two-dimensional case treated herein. 

Such a validation, based on the particular solution for the Kaplan 
section given in reference 7, has already been indicated in 
references 8 and 9. 

t f r hft method of analysis used herein is based on a conventional 
iteration procedure s imilar to that of the Bayleigh-Janzen method 
(reference 10) and the Acker et -Prandtl-Glauert method (reference 7). 
The Bayleigh-Janzen method expresses the potential as a power series 
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in the Mach number, whereas the Ackeret-Prandtl-Glauert method yields 
the potential as a power series in' a parameter characterizing the 
body thickness. The present method expresses the potential as a power 
series in a parameter that depends on all the physical parameters of 
the problem, namely, thickness ratio, free-stream Mach number, and 
specific heat ratio. The form of this dependence is not specified' 
initially. Instead, by going to the simultaneous limit of zero 
thickness ratio and free-stream Mach number unity in each step of 
the iteration procedure and retaining only highest-order terms, a 
power-series parameter (if it exists) is obtained as same combination 
of the previously mentioned physical parameters. The method thus 
constitutes a unification of the Rayleigh-Janzen and Ackeret-Prandtl- 
Glauert procedures in the small-perturbation transonic limit.. 

The preceding method of analysis turns out to be convenient for 
numerical application. Accordingly, as a by-product of the investi- 
gation, such application is made to the Kaplan section, and the 
results are compared with those of references 7 and 8. 'In particular, 
a value of the previously mentioned parameter of expansion of the 
power series for the potential is derived, which is estimated to repre- 
sent a boundary value between convergence and divergence of the series. 
The significance of such a boundary value is discussed in reference 11 
in connection with the so-called potential-limit phenomenon. 


AHALISIS 

General formulation. - The partial differential equation for 
the potential of an irrotational compressible, two-dimensional flow 
with free-stream velocity U is, in Cartesian coordinates x,y 
(fig. 1). 

[a 2 - (U+u^cp^ + (a 2 -^)^ - 2(Ufu)vq> X y = 0 (1) 

The following notation haB been used: (A more complete list of 

symbols is given in the appendix.) 

a local speed of sound 

U+u resultant velocity in x direction 

v resultant velocity in y direction 

cp perturbation potential, defined by u* cp x , v = qp y - 

The subscripts denote differentiation. 
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The local speed of sound a .is related to the free-stream speed of 
sound ag, the ratio of specifio heats 7, and the local velocity 
by the Bernoulli equation 


a 2 = ao 2 - ■2— (EUu+^+v 2 ) (2) 

In accordance with the Ackeret -Prandtl-Glauert type of procedure, 
equation (l) will be written in a form in which the linear terms appear 
on the left-hand side of the equation and the nonlinear terms on the 
right-hand side. A solution will be sought for high subsonic free- 
stream velocity and small lateral-distance ratio such that the flow 
pattern obtained by inclusion of the nonlinear terms will differ by 
only a small amount from that obtained with only the linear terms. 

The coefficient of qpyy in equation (l) is therefore expressed, 
with the aid of equation (2), in the form 



where 

Mq free-stream Mach number (=u/ag) > 

P 2 m 1 - Mq 2 (4a) 


^*M S Mo 2 + “o 2 ) 

For convenience the free-stream velocity is taken as the unit of 
velocity, so that u/tl and. v/CT may be written as u and v, 
respectively. The differential equation (l) can new be expressed 
in the form 


( P 2( Pxx' +c ^yy ) 

MoVV + ^(^x^j^xy 

The boundary conditions of the problem are that the perturba- 
tion velocities vanish at infinity and that the flow follow the 
contour of the body. Thus, at », 


1 - Mg 2 (2cp z +<p/^ y 2 ) = PM^qk+cftc 2 ) + ^ Mg 4 ?/ 


%cx + 


03 

t- 
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<P X - <Py. - 0 (6) 

The body is defined by 

7b ■ Tg(x) (7a) 

On the body, 

qp y a t(1+%)Sx (7b)- 


where 

t lateral distance ratio of body 

© 

g(x) function characterizing shape of body and together with 
its derivative g x is of order of magnitude unity 

and all lengths are expressed in teems of the chord o of the body 
as unit. 


The lateral distance ratio T, introduced in* equation (7a) to 
yield the order of magnitude unity for g(x), is a parameter to be 
regarded as including both thickness ratio and angle of attack. The 
parameter t reduces to thickness ratio for a symmetrical airfoil 
at zero lift, to angle of attack for a flat plate with lift, and to 
camber ratio for a fore and aft symmetrical curved plate at zero 
angle of attack. 

In order to obtain the Laplaoian on the left- h a nd side of 
equation (5), the affine transformation 

11 = Py (8a) 


F(x,n) = -§ <P(*,y) 

r 


(8b) 


is introduced. The transformation from qp to F in equation (8b) 
is made in order to incorporate the factor F y/P 2 indicated by 
equations (5) and (8a) into the solution for the potential. The 
differential equation (5) thus becomes 
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" m W p2 V) " + 

(®* + fj '** ^ Mo 4 ^ 2 ) JW« 0 2 r 4 V*W, + 

*£'*)**» (9) 

where AF Is the Laplaolan 

\ 

AF =F XX + (10) 

3ii terms of the new variables, the "boundary conditions (6) and (7) 
"become at 


F x = = 0 (11) 

and on the "body (subscript b) 

% - 0 T g(x) (12a) 


T gr 

^ - nr F * 




(12b) 


The formulation of the problem is thus far exact. A solution 
is now sought applicable in the limiting case T— >0, 0 — >0j this 

limit will be referred to as "the small-perturbation transonic limit,” 
* In the range of values of T and p considered, the perturbation 
velocities will be assumed small compared to free-stream velocity, 
or |u| , |vj«l. The right-hand side of equation (9) is thus con- 
sidered as producing a small perturbation from the linear case so 
that a solution to the systems (9), (11), and (12) will be sought 
in the form 


o' 


12 3 

3?=3?+3? + 3?+« * . (13) 

in which each term is of a lesser order of magnitude than the 
preceding one. Accordingly, as the equivalent of boundary condi- 
tions (11) and (12) the following boundary conditions on 
1 2 

P, V , • . will be taken at » : 
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II n 

*ac ■ ° °> n = 1, 2, 3 . . . (14) 

on the "body 

I!* = 3 Tg( 2 i) (15a) 

^ t ^ Tr™ 

“ p Sx “ "^ 3 “ Sx (15b) 


n T n 

F il “ p Sx F x = 0 Ca j* 1) (15o) 

In order to obtain a solution of the systems ( 9 ), (14), 
and (15), equation (13) is. inserted into the differential equa- 
tion (9) and a typical Laplacian term, on the left- such as 
n‘ 

AF, is equated to the sum of those terms on the right that 
satisfy the following two conditions: (a) The superscript n-1 is 

present; and (b) any or none of the superscripts n~ 2 , n -3 ... 1 
are present. The resulting solution of the non-hamogeneous Laplace 

equation for F consists of a sum of terms of which, for a range 

of Mq near unity and T near zero, same will be of highest order 

o n 

of magnitude. These terms constitute the solution for F in the 
small-perturbation transonic limit. It will be shown that if T 

and • P are allowed to approach zero in a particular manner then 
, n 

the remaining terms will vanish and the solution for F will have 
a unique finite limit. 


Transformation to elliptio coordinates. - The Cartesian coor- 
dinates z,tj are suitable for obtaining the flow past a two- 
dimensional profile that is periodic with respect to x (for 
example, a wavy wall) but are inconvenient in problems involving 
Isolated airfoils. In order to obtain explicit solutions fear the 
flow past an isolated two-dimensional body,, it is convenient to 
introduce a system of elliptic coordinates (reference 12 , p. 156). 
The transformation from Cartesian coordinates x,rj to elliptic 
coordinates s,t is given by (fig. 1 ) 

x + Itj = cosh (s+it) (16) 

The real and imaginary parts of equation (16) yield 
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x = cobIl s cos b (17a) 

tj = sinh s sin t (17b) 

This curves s = constant, t = constant are confocal ellipses and 
hyperbolas, respectively, with, the common foci at x = ± 1, q = 0. 
The values of s range from. 0 to o=, whereas t varies between 0 
and 2it. The Laplacian AF is, in elliptic coordinates. 


AF 


■®ss + -^tt 
J(s,t) 


O 

(18) 


where J(s,t), the Jacobian of the transformations (17), is given 
by 


J(s,t) = XsTi t - x^ 

= sinh? s + sin 2 t 
= cosh 2 s - cos 2 t 


(19a) 

(19b) 

(19c) 


Transformations from derivatives with respect to x and tj 
to derivatives with respect to s and t will be needed for the 
subsequent analysis. From the transformation equations (17) these 
relations are 


F x = (sinh s cos tf s - cosh s sin t F t ) (20) 

(cosh s sin t F s + sinh s cos t Fj.) (2l) 


First approximation. - The solution for the first approxima- 
tion F satisfies the homogeneous Laplace equation 

1 1 1 

AF = P BS + F tt = 0 (22) 

anfl the boundary conditions (14), (15a), sued (15b), which in 
elliptic coordinates are at <» 
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1 . 1 1 . 

jjr~y (sinh s cos tF s ■ cosh s sin t F^) = 0 

1 , 1 1 

(cosh s sin t F s + sinh s cos t F t ) = 0 

on the body 


% = 3Tg(2:) 


(23a) 

(23b) 


(24a) 


1 , ^ 1 1 

(cosh s sin t Fg + sinh s cos t F t ) - 

f ^ jfa/fc) ( sinh s 008 * F a ~ OOBh s sin t F t ) = lEM g 

p 3 * 

(24b) 

A solution of equation (22) that satisfies the boundary condi- 
tions (23) is 


1 

F 



1 

(Ajj sin nt 


a 

+ Bg, cos nt)e“ ns 


1 

+ct 


(25) 


It is necessary to exclude the positive exponentials e 110 from, 
equation (25) because, from equations (23) and (19b), the order of 
11 2 
F Q and F t at Infinity must be less than e s . The constants A n 

and B n in equation (25) are determined from the boundary condition 

on the body (equations (24a) and (24b)). The constant C is sub- 
sequently determined to yield a circulation around the body by a 
Kntta-type condition (such as rear stagnation point at trailing 
edge). Small-perturbation analysis usually neglects the term in 

F x occurring In equation (24b). This term will therefore be assumed 
negligible, an assumption to be Justified subsequently on the basis 
of the resulting solution. 

In obtaining the small-perturbation transonic limiting solu- 
tion for F, equation (25) is inserted into the limiting form, for 
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&Ij — > 0, of equation (24b) with, the F x term neglected. The 
limiting forms of equations (19), (20), and (21) as s — »0 are 


11m J(s,t) «= sin 2 t 
s — *0 


11m F x = - 
s — >0 


*t 

sin t 


F, 


llm F- = — — : 
. - *i sin t 

s— j0 


and the resulting equation for determining the constants A_ 
1 

and Bn 1 b 


(26) 

(27) 

(28) 



1 1 tIm 

n^ sin nt + cos nt) g- ^(x) sin t 


(29) 


The more rigorous alternative process of Inserting equation (25) 

into equation (24b) with the F x term dropped and only then going 
to the limit — >0 yields the same equation (29) . Evaluation of 

the boundary condition for F on the axis, that is, by equation (29), 

1 

will therefore yield the correct limiting solution for F. It may be 

noted that in order to obtain a solution for F it is necessary that 
the quantity y/ P 3 la equation (29) be finite in the small- 

perturbation transonic limit t— >0, 6 — *0. 


Equation (29) does not contain the variable s inasmuch as 
the limiting form of g^x) as s — *0 is understood to be sub- 
stituted Into equation (29). The right-hand side of equation (29) 
may now be expanded Into a Fourier series in t and the coef- 
ficients of corresponding functions on both sides equated. The 

1.1 

resulting solution for A n and may be indicated as 
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(30a) 

1 1 

Bn = Vb)K 

(30b) 


where 


t£m 

P 3 


(31) 


and a^b), b n (b) are functions of the "body shape. (The symbol b 
within parentheses denotes the numerical coefficients characterizing 
the body shape, including angle of attack, that arise in expressing 
gg(x) as a trigonometric series in t.) numerical evaluations 

such as for the Kaplan section to be made herein Indicate that 
11 

%0>>> b n (b) are of order of magnitude unity. 

The solution for F is thus by equations (25) and (30) 


1 

F = K 


1 1 na 1 
(ajj sin nt + b n cos nt)e _ns + ct 


(32) 


1 i 1 

in which the substitution C = Kc(b) ^has been made, where c(b) 

may be considered of order unity if a^b), b^b) are of order 
unity. 

1 1 

The functions and bjj do not depend on the coordinates s-j, 
because they were obtained from equations evaluated on the axis of 
the body (s b = 0). The reduoed potential function F therefore 
depends an s only through the exponential e _Ils , which is nearly 


s 


0. 


constant^for small values of s and approaches unity for 

Hence, F approaches a finite limiting solution as T and p both 
approach zero in such a way that K is kept finite. This limiting 

solution for F holds for a class of profiles with a common- shape 
and' for a fixed value of K, so that a similarity 


factor g(x) 
rule exists. 


If the limiting solutions for T, p— >0 in the 
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higher approximations can likewise he expressed in terms of the one 
parameter K, it will then he possible to say that a transonic 
similarity law exists for two-dimensional flows. 


The horizontal and vertical disturbance velocities are, from 
equations (8), (20), (21), and (32), 


1 q2 1 T 1 — — 

<Px = f ffjc* 5 - |f jfs/fc) ^ sinl1 3 003 t ^ / n ( a n sin nt+h^ 008 nt)e _DB + 


cosh, s sin t 


n(a a cos nt-h^ sin nt)e _ns + c 


(33) 


1 a 3 1 T 

9 y =^- “ j^eT^tT ^ oosl1 0 sin 11 ^ sin nt+bjj cos nt)e - ns - 


M 


£n(i 


sinh s cos t 


n=0 


1 1 1 
n^ cos nt-bjj. sin nt)e”“ B + c 


(34) 


The total disturbance velocity A is defined by 


1 

A = 


(1-kd) 2 +q> 


- 1 


(35) 


The trigonometric and exponential functions occurring within the 
summation signs in equations - (33) and (34) are of cutler unity and 
it is reasonable to assume that the results of summation will be 

1 1 

of the same order. The ratio of cp_ to <p_ is of order P so 
' 1 2 7 1 
that the quantity <p v is negligible with respect to cp z in the 

7 1 1 
small-perturbation case cp z «1. The disturbance velocity A is 

therefore approximated by cp z . In the first approximation given 
by equations (8), (31), (33), and (34), the similarity rule is 
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therefore simply the veil -known proportionality of the resultant 
perturbation -velocity to thickness ratio and angle of attack 
(represented by the symbol t) and to l/3, vhich is the 
Prandtl-Glauert rule. 

2 

Second approximation. - - The second approximation F is a 
solution of the non-homogeneous Laplace equation 



which satisfies the boundary conditions (14) , (15a), and (15c). As 
previously mentioned, in solving equation (36) only those terms on 
the right-hand side will be retained that are of highest older of 
magnitude in the small-perturbation transonic limit T,p— >o, K 
finite. It is thus first noted that the term involving the 
/I 1 \ 

Laplaoian factor I F sa +F .^ J vanishes by virtue of the solution 

for the first approximation. Seoondly, equations (20), (21), 

1 

and (32) show that F and its derivatives are of order of magni- 
tude K. In the numerical applications hereinafter given, the 
permissible range of I, is approximately 0<K<0.5. It is 
therefore convenient to assign the order of magnitude unity to K. 

The order of magnitude, as 3 — >0, of any term on the right-hand 
side of equation (36) therefore depends on the power of 3 occurring 

1 1 

in that term, and only one term, involving F z F^, does not have a 

factor of 3 to some positive power. In the small -perturbation 
transonic limit, all terms on the right-hand side of equation (36) 

1-1 . . 
except the term in F_ F^ are therefore negligible. Equation (36) 

2 

for the second approximation F thus simplifies to 


2 2 
F bb + P tt 


2j(s,t)F 3 F Z2 


(37) 
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For the solution of equation (37), F i S first expressed as 


F = KE(s,t) 


where, "by equation (32), 


B(s,t) = (a^ sin nt + bn cos nt)e _ns + ct (38) 


The x derivatives of E are given by, from equation (20), 
i (sinh s cos t Eg - cosh s sin t Ejj) 


(39) 



s cos 2 tR ss + cosh 2 s sin 2 tRj-j. - 


2 cosh s sinh s cos t sin tSg^. + 


sinh s cosh s 


1 - £ . °P* . * (Btth? B - Bln 2 t) Eg + 


cos t sin t 


1 + ^ — c . °gj 1 . - 8 . (sinh 2 s - sin 2 t) | BfcJ* (40) 


In terms of B, equation (37) becomes 


F ss + F tt - HK 


(41) 


where E^ and B^ are obtained from equations (38) to (40). The 
right-hand side of equation (41) is now to be expanded in a Fourier 
series in t. A typical term of equation (41) has a singularity 
at the end points t = 0,st for s = 0, because of the factor l/j, 
which becomes Infinite at these points. Because the limiting solu- 
tion for F for s — >0 is of principal interest, it may be neces- 
sary to exclude a small and finite^region around the end points 

from the domain of def ini tion of F. This objection is not 
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serious, however, because the end points must in any event he 
excluded from consideration if stagnation points, which violate 
the assumption of small perturbations, occur there. When the 
right-hand side is expanded in a Fourier series in t, equa- 
tion. (41) becomes 


bj^s) sin nt + E^(s) cos ntj (42) 

where h^s) and h^s) are functions of s. 

2 

A solution far F is now sought in the form 


ss 


+ *tt = e2 


2 

F 


03 



s) sin nt + q^s) 


cos nt 


■] 


(43) 


where q n (s) and q^s) are functions of s to be determined. 
Insertion of equation (43) in equation (42) yields 


a> 

j(q n "-n 2 q n ) sin nt + (q n ,, -n 2 q n ) cos nt] = ^ (1^ sin nt + ^ cos nt) 

(44) 


where the primes denote differentiation with respect to s. 
Because equation (44) must hold for each value of n and for all 
values of t, the coefficients of corresponding trigonometric 
functions on both sides may be equated, which yields the ordinary 
differential equations 


In” ~ n2 9n * h. ( 45a ) 

q n ” - n 2 ^ = h^ (45b) 

Equations (45a) and (45b) are of the same type so that attention 
may be confined, say, to equation (45a). The solution for q n (s) 

2 

consists of a complementary function, taken as a^e to satisfy 
the boundary conditions at infinity, and of a particular -integral 
th at may be expressed in terms of two indefinite integrals by 
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r n( s ) 



hn(s)ds 


(46) 


In order to determine the behavior of the particular integral 
r n (s) as s— *0, the form, of hj^s) is first examined. Because 
a region around the end points is excluded if stagnation points 


occur, the quantity 


g 

cosh s 


occurring in equation (40) may be 


expanded in an absolutely convergent power series in 


cos 2 t 

g < 

cosh s 


The 


right-hand side of equation (41) may now be expressed as a Fourier 
series in t. The function h n (s) may therefore be expressed as a 


sum of terms each of which has the form cosh? 5 - se _ms or 

sinh s cosh k se _ms where k and m are integers, m positive, and 
each term is mult i plied by a coefficient that depends an a body 
shape and is of order unity. The particular integrals corresponding 
to these forms are, from equation (46), expressible as functions 
also having these forms. It then follows that the function r n (s) , 
corresponding to hu(s), and its derivatives r n *(s), ^"(s) are 
of order of magnitude unity and have a finite limit as s — ^0. 


Writing r n (s) as the particular integral corresponding to b^s) 

u 


gives the following solution for F: 


1 =^ 


A circulation term ct Is included in equation (47) in order to 
satisfy the Kntta condition as in the first approximation. The 
2 2 

coefficients a^ and bn are determined from the boundary condi- 
tion on the body (15c) which, with neglect of the F x term and by 
equation (21), is in elliptic coordinates 

2 2 

cosh s sin t F a + sinh s cos t F^ = 0 (48) 


2 2 
(a^e^+rj!) sin nt + ? (bne^+ra) cos nt + ct 


(47) 


The result of combining equations (47) and (48) and going to 
the limit s^ — >0 is the same equation as that obtained from 
combining equation (47)’ with the limiting form, as — * 0, of 
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equation (48). It Is therefore permissible, as in the first 
approximation, to evaluate the boundary condition on the axis of 
the body for vhich equation (48) becomes 


P s = 0 


( 49 ) 


Substitution of equation (47) in equation (48) yields 


OB OO 

2 V'' 2 

(-najj+rjj*) sin nt + 2 , (“nb n +r n ’) cos nt 
n=0 n=0 


= 0 


(50) 


For equation (50) to hold for each value of n and all values 
of t, it is necessary that 


z 

-n&ja + r n * «= 0 



(51a) 

2 

-rib„ + r * = 0 
n n 



(51b) 

Solving equations (51a) and (51b) for the coefficients 
yields 

2 

2 

and b n 

ii 

1? 



(52a) 

2 . 
*n c r n'/ n 



(52b) 

Equations (52) show that the coefficients 
order of magnitude unity. 

2 2 
an and 

are 

of. 


The small-perturbation transonic limiting solution for the 
* 2 ? 
second-approximation potential F is thus proportional to K 
and to a function of body shape and position of order of magnitude 
unity, which approaches a finite limit as s— >0. The derivatives 
2 2 

of F and, in particular F z , also have these properties. From 
the results obtained thus far it is clear that similar solutions 
will be obtained in succeeding approximations because the right- 
hand side of equation (9), in the small -perturbation transonic limit, 
will always, in any approximation, be made up of power of E and of 
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functions of order of magnitude unity that approach a finite limit 
for s — >0. The various approximations may he conveniently arranged 
in the form of a power series in E and the. final results indicated 
as follows : 

On the body s = s^ ~ 0 


= -*?j A = a^K + SgE^ + a gE 3 + • » « 
In the flow field s ~ 1 

F* ■ v = b x E + h 2 K 2 + bjE 3 + . . . 


(53) 


(54) 


where 

A disturbance velocity an body 

v disturbance velocity in flow field 

a^ag ..., b-j^bg ... functions of body shape and position of order 

of magnitude unity 

The velocity on the body and in the flow field in two-dimensional 
transonic flow past a thin airfoil therefore depends on Mach num- 
ber Mqj specific-heat ratio 7, and lateral distance ratio T 

only through the combination K = The small-perturbation 

transonic limiting form of the parameter E is iT/p 3 , where 

r s Tim r = This limiting form of the parameter E is 

%- >1 M 2 

a simple function of the similarity parameter of reference 1. Thus 
the transonic similarity law derived here agrees, in the small- 
perturbation transonic limit, with that of reference 1. It is 
apparent that the validity of the derivation of this law depends 
principally upon the properties of the potential in the neighborhood 
of the axis. 

F inally , the previous assumption that in the body boundary 
condition (12) the F x term could be neglected with respect to 
the F-jj term will be Justified. Because F x and Fq are of the 
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same order, the ratio of the F z term to the team in equa- 
tion ( 12 b) ie of order t/P. For a transonic flow with the 
parameter K of order one, t is of order P^. ratio t/P 
is therefore of order P 2 , which is small compared to unity in 
the transonic limit' p — >0. Bence, the neglect of the F z term 
was Initially permissible. 

ILUJSTEATI7E EXAMPLE 

Although the prime purpose of the method of analysis was to 
verify the transonic similarity law, the method of analysis may 
now be used to calculate flows past specific airfoil sections. As 
an illustration, the first two approximations for the flow past a 
Kaplan section will be obtained in the limiting form given by 
equation (53) . The parametric equations of the Kaplan section are 
(reference 7) 


x-jj = cos 9 


■£ (cos 0 - cos 39) 


y-fa = ^ (3 sin 0 - sin 30) 

where the parameter 0 varies between 0 and 2x and the coordinates 
x b, y b, 111 conformity with the convention of reference 7, refer to 
the semichord as the unit of length. 

First approximation. - The slope dy^/dx^ is given to the 
first order in t by 

^ y b 37 / cos 9 - cos 39 \ 

dx-jj 4 \ sin 0 ) 


= -3t sin 0 cos 0 

In terms of x^ the slope is then given, to the first order in t, 
by 
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° In general, the slope Is, from equation (7a), 


fob 

^b 


- T ^(X) 


so that here 

83 .( 2 ) = -3^ ^l-x b 2 (56) 

The coordinates x^/lb are given In terms of elliptic coordinates 
by equation (17) 


3 b » cosh s cos t (17a) 

T)b = sinh s Bin t (17b) 


Because the boundary condition is to be evaluated at s 
tian (17) becomes 


Xfc as cos 
^b c 0 


t 


1 


0 , equa- 


(57) 


The function g^x) is therefore given by, at s = 0, 


gj.(x) = -3 cos t sin t 


(58) 


Equation (58) could have been obtained directly from the 
expression for g x (x) in terms of 6 by noting that, in the 
limit s — > 0 , the parameter 9 becomes identical to the 
coordinate t. 


1 

The first approximation F 


op 



1 

(An sin 


is, from equation (25), 
nt + Bjj cos ntje ” 118 


(25) 


The boundary condition for F is, from equation (29), 
00 

sin nt + Bjj oos nt) = -Eg^x) sin t 


n =0 


(29) 
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parting the -value of gj.(x) tram, equation (58) in equation (29) 
yields % ' 


OO 

1 1 

/ , n (^n sia n-fc + ^n cos n-fc ) = — (cos t - cos 3t) 
n =0 4 

11 

The coefficients and Bq are therefore given hy 


1 

A n c 0 


% -f X 

1 T 

B 3 = - i K 

ki - 0 , n { 1 , 3 J 


(59) 


The first approximation F is t.Vm a 


I s « j (3e" s cos t - e“ 3s cos 3t) 


The limiting form of F for s — >0 is 


(60) 


1 ^ 

F o — (3 oos t - cos 3t) . (Si) 

The limiting form of the velocity A for s -*0 is, hy equa- 
tions (27), (33), and (35), ' 


1 

C M " " f K 008 2t 


(62) 


The m a xi mum value of A occurs at the 'midchord (t = «/ 2 ) and is 
given hy 


r iAlflx/f 2 =f K 


(63) 
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2 

Second approximation. - The second approximation F is the 
eolation of the nan-homogeneous Laplace equation ( 37 ) 

*ss + *fct = 2J (BjtfrJzx ' (37) 

1 

The x-derivatives of F may he obtained hy successive applications 
of equation (20). It is more convenient, however, to first intro- 
duce the conjugate complex variables 


£ = s + it 

► 

£ » s - it 

Equation (60) becomes, in terms of £ and T, 


(64) 


F « | ^3e“£ - e“ 3 5 + 3e“£ - e" 3 ^ (65) 

The derivatives with respect to x are given in terms of deriva- 
tives with respect to £ and £ by 


F 


F = 
■ C X 


sirih 


£ sinh £ 


( 66 ) 


■2X 


1 F>>- S2fei- F. + l__F 7r - . c -? s . h -.^_ Fr + 

sinh 2 £ ^ sinh 3 £ ^ sinh 2 £ ^ sinh 3 £ * 

2 


F-, 

sinh £ sinh £ * * 

The x-derivatives of F are, from equations (65) to (67) 


(67) 


i = -| E (e- 2 S + e^) 


*xx “ § * ' 9 


sinh £ sinh 


0 

1 

N, 
<r»l 
' — -• 

(68) 


(69) 
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The Jacobian J(£,£) and the Laplacian AF(£, £) are given in 
terms of £ and f iy 

J(£>£) - sinh £ sinh J (70) 

AF(£,I) = i il-r (71) 

sinh £ sinh. J ^ 

, » 2 
The differential equation (37) for the second approximation F 

becomes 

2 _ 

4F^ = - | K 2 (e" 2 £ + e" 2 £)(-2£ sinh £e + -2$ sinh £e) (72) 

If equation (72) is transformed hack to s,t coordinates, there 
is obtained 


F sb + <= ^ K 2 J^(e“ 5s - e“ 3s ) cos t + e“ 3s cos 3t - e“ 3s cos 5tJ 

(73) 

In order to solve equation (73), a solution is assumed of the 

form 


2 9 9 T 2 2 2 "I 

I o j F <Li(s) cos t + q 3 (s). cos 3t + q 5 (s) cos 5t (74) 

Insertion of equation (74) into equation (73) yields the following 
ordinary differential equations for the q functions: 


*1 - *1 


e 


-5s _ 3 -3s 


2 „ 2 ,- a 
a 3 - 99.3 - e" 5& 

l 5 ” - 2s| 5 = -e" 3s 


The solutions for the differential equations (75) are 


( 75 ) 



Kr)CO 
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2 , , 2 _ s • 2 , ■ 
^(s) = A^e B + r-j^s) 

2 2 2 
0.3(3) = A 3 e- 3s + r 3 (s) 

tgte) = A 5 e- 5s + r 5 (s) 




(76) 


2 2 2 

The particular Integrals r^s), r 3 (s) and r 5 (s) may he evaluated 
hy means of equation (46) 


r n (s) 


= e” ns J ] 3 2118 ds /V“ h^(s) ds 


d 


(46) 


The result is 


*i(») - 3 e " Ss - | 8 ' 38 
*3 (s) ■ i °" 5s 


*s< B * “ TO e ' 3B 


(77) 


2 2 2 

The constants A^, A 3 , and A g are determined from the boundary 
condition (49) which, when applied to equation (76), yields 


2 . 

A x = 1/6 


Ag = -5/48 f 

' * 

'5 


2 , / 

A_ = -3/80 


(78) 


The solution for F is therefore 


| e2 [(I e_a -f e ‘ Ss+ 3 e ' SB ) C0B 4 + (S e ‘ 3s + TO e ' S °) 00B 34 + 

(SO °' 5S + M ^ °° S 5t ] (79) ’ 
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The limiting form of F as s — *0 is 


2 / \ 

F = x K 2 f 008 'fc - — =r cos 3t + JL cos 5t ) 
4 \12 24 40 / 


( 80 ) 


The limiting form of the velocity A as s— 40 is, "by equa- 
tions (27), (33), and (35), 


' 2 

2 r M A 9 K 2 /l l i \ 

4 HrTt (n -t " - sia 3t + i sin 5t] 


Fv = 




8 


8 


(81) 


The maximum value of A occurs at the midchord (t = jc/2) and 
is given by 


^WWx D 3 k 2 


4 


(82) 


CQMPAEISOn OF HESUIffiS WITH KAPLAH 

The first and second approximations for the velocity potential 
for the Kaplan section as obtained by the Ackeret -Prandtl -Glauert 
iteration method are given by equations (B-17) and (C-5), respec- 
tively, of reference 7. The limiting forms of these expressions 
as T — >0 and P—4 0 are in agreement with the present results, 
(equations (60) and (79), respectively). Hence the legitimacy in 
the present report of obtaining a limiting solution by going to 
the limit in the differential equation rather than in the final 
solution is verified. The first procedure requires, of course, 
much less work than the second. It may now also be inferred that 
the limiting forms of Kaplan's results in the third approximation 
are the same as would be obtained by the present method in the 
third approximation. In particular, the third^approximation for 

the maximum velocity increment (at mid chord) A mB7 , given by equa- 
tion (17) of reference 7, has as limiting form 

3 

^WWx 87 ^3 

pr = 80 E 


( 83 ) 



26 


NACA IB 2191 


Tha small-perturbation transonic limiting form for the ma-Timum 
velocity on the surface of the Kaplan section is therefore given, 
to three approximations, by 



3 

4 


r ,87 e 2 
K + 80 K 


(84) 


The maximum velocity increments given in similarity form by 
equation (84) an d the corresponding values obtained from refer- 
ence 7 for several thickness ratios and Mach numbers are plotted 
in figure 2 as a function of K. It is seen that the results for 
the finite thickness ratios do not differ negligibly from those for 
the limiting solution, so that same error would be incurred in 
extrapolating the limiting solution for zero thickness ratio to a 
finite thickness ratio. For airfoil sections mare closely approxi- 
mating the ellipse in shape, however, the effect of thickness ratio 
on velocity may be expected to be smaller than for the Kaplan 
section. This conclusion is indicated in unpublished calculations 
by the method of reference 11 and mentioned in presenting 
reference 9. 


The local Mach number M corresponding to a given velocity 
increment may be obtained from Bernoulli’s equation which, in 
small-perturbation form, is given by (reference 11) 


M 2 -l = -P 2 + 2r M A (85) 

Equation (85) can be expressed in the similarity form. 

Curves of constant value of the’ parameter (M 2 -l)/p 2 in equa- 
tion (86) have been plotted in figure 2 to indicate the local Mach 
number corresponding to a given velocity increment. 

It may be noted that the small-perturbation transonic limiting 
value of K derived in reference 8, corresponding to the first 
appearance of locally sonic velocity on the Kaplan section, is given 
by the intersection in figure 2 of the t = 0 curve with the 

(M 2 -l)/p 2 = 0 local Mach number contour. In order to obtain local 
Mach numbers corresponding to the finite-thickness-ratio velocities 
in figure 2, the exact form of Bernoulli’s equation should be 
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used, rather than the small-perturbation form that leads to equa- 
tion (86). Accordingly, contours of local Mach number unity, 
corresponding to A m , based oh the exact form of Bernoulli’s 
aquation are shown in figure 2(b) for several thickness ratios. 

In order to aid in obtaining sonic values of E for other thick- 
ness ratios, the interpolation curve AB has been drawn through 
the intersection of these contours with the corresponding velocity 
curves. Local Mach numbers other than unity for finite thickness 
ratio may be estimated by comparison with the various contours of 
local Mach number unity. 


Each velocity curve of figure 2(b) has been terminated at a 
particular upper value, of K, indicated by a circle. These 
values of K, called Ej, were obtained from the values of M 0 j 
for given t derived in reference 7. The value of Mq j is the 
limiting Mach number above which, for the given T, the series 
solution for the velocity may be expected to diverge. The value 
of Ej for the limiting solution T — *0 was estimated from equa- 


tion (84), as in reference 7, by comparison with the corresponding 


terms of the harmonic series 



The ratio of the terms of 


the series solution to the corresponding terms of the harmonic 
series indicated a value of K below which this ratio formed a 
decreasing sequence. For a discussion of the significance of these 
so-called potential limit quantities, see reference 11. 


For convenience in converting the similarity parameters in 
figure 2 to more directly usable quantities, such as velocity 
increment A, a plot of the similarity parameter K against 
free -stream Mach number Mq for several values of the thickness 
ratio T is presented in figure 3. 

The aotual values of the velocity increment A, free-stream 
Mach number Mq, and looal Mach number M at the potential limit 
anfl the first two of these quantities at the critical Mach number 
Mq or (local Mach number unity) obtained from the data of refer- 
ence 7 and shown in figure 2 are presented as functions of thick- 
ness ratio T in figure 4. For comparison are included the 
corresponding values derived from the similarity parameters fear 
T = p = 0 given also in figure 2. The local Mach number for the 
l imiting solution should, for consistency, be calculated from the 
small-perturbation form of Bernoulli’s equation (85). In order to 
Indicate the approximation introduced by equation (85)-, the local 
value of M for both the critical and potential limit conditions 
has been also obtained from the exact form of Bernoulli’s equation. 
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In order to compare the exact and limiting solutions with: regard <» 

to flow-field conditions, calculations have been made of the maximum to 

lateral extent of the local supersonic region at the potential limit 
Mach number for several thickness ratios by means, of both solutions. 

The results for the exact solution were obtained by the method of 
reference 7 (equation (18)), whereas those for the limiting solution 
were obtained in 'a similar manner by equations (8a), (60), and (79). 

As in reference 7, only the first two approximations were used in 
making the calculations, because, as noted in reference 7, the 
disturbance velocity decreases rapidly with distance from the body. 

The results are shown in figure 5, in which the maximum height in 
chord lengths of the supersonic region has been plotted as a function 
of thickness ratio. The results for the limiting solution are in 
surprisingly good agreement with those for the exact solution, 
especially in view of the rapid change of lateral extent of the 
local supersonic region with Mach number in the neighborhood of the 
potential-limit Mach number shown in references 7 and 13. It is 
noted that the height increases rapidly with decrease in thickness 
ratio at the small thickness ratios. This behavior is in conformity 
with the similarity law, which gives a fixed value of the nondimen- 
sional parameter corresponding to the height of the supersonic 

region. The parameter Pj approaches zero as the thickness ratio 
goes to zero so that the height y^ = q^/pj (equation (8a)) cor- 
respondingly becomes Infinite. 

The results given by Kaplan far the circular arc (reference 14) 
may be analyzed in a similar manner. The maximum velocity increment 
obtained from equation (18) of reference 14 yields, in the small- 
perturbation transonic limit. 


— = 4 + 6K + M E 2 (87) 

h 3 

where h is the camber ratio (ratio of maximum height to chord 
length), and K = hT M /p 3 . The maximum velocity increments given 
by equation (87) and reference 14 are plotted in similarity form 
in figure 6. The velocity increment A and free-stream Mach 
number Mq at the critical Mach number Mq qj. for both the exact 
and limiting solutions are presented as functions of h in figure 7. 
In the same manner as for the Kaplan section, the local Mach number 
obtained from the exact form of Bernoulli’s equation is included. 

The values of for the exact and limiting solutions are in 

better agreement for the circular arc than far the Kaplan section. 
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coucnroniG eemaeks 

A transonic ^similarity law for two-dimensional flow has "been 
derived by von Harman; in the derivation of this law, the assump- 
tion has been made that the boundary condition at the body may be 
satisfied on the axis. This assumption has been investigated herein 
by an iteration procedure .similar to the Bayleigh-Janzen and 
Ackeret -Prandtl-Glauert procedures, in each step of which the boundary 
condition is satisfied exactly an the body. The resulting solution 
is of the form, required by the existing transonic similarity theory, 
thereby validating the boundary-condition assumption. The question 
of the effect of stagnation points upon the similarity law, however, 
has not been Investigated. 

In using the method, the assumption was made that the small- 
perturbation transonic limiting form of a solution of the governing 
differential equation is the same as the solution of the limiting 
form of the differential equation. This assumption was verified by 
comparison with the exact solution for the Kaplan section. 

The present method of analysis offers seme advantages in 
simplicity and amount of computational work as compared to the exact 
method in applications to specific airfoils. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland > Ohio, May 24, 1950. 
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APPENDIX - SYMBOLS 


The following symbols are used In this report: 
Aj^BqjC arbitrary constants (equation (25)) 
a local speed of sound 

a^bj^c function of body shape (equations (30) and (32)) 
&q speed of sound in free air 

c chord of body 


P 

g(x) 

h 

^n^ 

J 

K 

Mq 

E 


transformed velocity potential 



function characterizing. shape of body 
camber ratio 

functions of s (equation (42)) 
Jacobian 


similarity parameter 



free-stream Mach number (= U/a Q ) 
functions of s (equation (43)) 
function defined by equation (38) 


r n particular integral (equation (46)) 

s,t elliptic coordinates 

U free-stream velocity 

u disturbance velocity in x direction 

v disturbance velocity in y direction 

x,y Cartesian coordinates 

p fljl - Mq 2 


1379 



HACA 03T 2191 


r 2+i 

2 

r M Mo 2 + Mo 2 ) 

7 ratio of specific heats 

A Laplacian 

J (= b + it) 

? (- B - it) 

Tj transformed y- coordinate 

A disturbance velocity on body 

v disturbance velocity in flow field 

T lateral distance or thickness ratio 

tp disturbance velocity potential 

Subscripts : 
b on the body 

or critical 

Z potential limit 

max maximum 

The symbols x, y, tj, s, and t used as subscripts indicate 
differentiation with respect to that variable. 
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Similarity parameter, K 
(a) First and. second approximations. 

Figure 2. - Variation of maximum velocity Increment vith similarity parameter 

for Kaplan section* 



Maximum reflooity increment 
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Similarity parameter, K 

(h) First three approximations. 

■ » 

Figure 2. - Concluded. Variation of maximum velocity Increment vith similarity 

parameter for Kaplan section. 
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Figure 4. - Variation of maximum velocity increment, free-stream Mach number, 
and local Mach number vlth thickness ratio at critical Mach number and poten- 
tial limit for Kaplan section. 


1379 



Maximum height of supersonic region 


MCA TS 2191 



Thickness ratio, T 

Figure 5. - Variation of Tnarimnni height of supersonic region at potential limit 
vlth thickness ratio for Kaplan section. 
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Figure 6. - Variation of maxim u m velocity increment vith. similarity parameter 

for circular arc section. 
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Figure 6. - Concluded. Variation of m a xi m u m velocity increment vith s imilar ity 

parameter for circular arc section. 
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